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Abstract 

We study the change in condensation energy from a single pair of fermionic atoms to a large number 
of pairs interacting via the reduced BCS potential. We find that the energy-saving due to correlations 
decreases when the pair number increases because the number of empty states available for pairing gets 
smaller ("moth-eaten effect"). However, this decrease dominates the 3D kinetic energy increase of the 
same amount of noninteracting atoms only when the pair number is a sizeable fraction of the number of 
states available for pairing. As a result, in BEC-BCS crossover of 3D systems, the condensation energy 
per pair first increases and then decreases with pair number while in 2D, it always is controlled by the 
"moth-eaten effect" and thus simply decreases. 
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1. Introduction 

It was known for a long time that a 3D sys- 
tem with a weak attractive potential cannot sus- 
tain a bound state. Cooper however showed that 
in the presence of a frozen Fermi core, a pair of 
electrons with opposite spins can form a bound 
state with zero total momentum, no matter how 
weak the attraction is.l.]. Note that this "sin- 
gle pair" state already is a many-body state be- 
cause, even if the frozen core electrons do not in- 
teract, they still, because of Pauli blocking, pro- 
vide a finite density of states which is of im- 
portance for pairing. Turning to more than one 
pair is difficult due to the Pauli exclusion prin- 
ciple between a fixed number of paired electrons. 
To overcome this difficulty, Bardeen, Cooper and 
Schrieffer proposed an ansatz for the many-body 
state in the grand canonical ensemble - with pair 
number not fixed - which in the presence of a 
frozen core, leads to an energy lower than the 
free electron energy, even in the limit of an ar- 
bitrarily small potential^. Gor'kov and Melik- 
Barkhudarov then showed that the frozen core is 
not mandatory, provided that one uses a renor- 
malized attraction measured through the low en- 
ergy scattering amplitude [1]. Later on, EaglesQ, 
LeggettQ and also Nozieres Schmitt-RinkQ ex- 
tended the BCS idea to bridge molecular BEC 
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with Cooper pairing. To do it, they vary the 
potential amplitude while using a BCS-like grand 
canonical wave function without frozen core, i.e., 
all k states are involved. Their work raises the 
complementary question: how, when the potential 
is fixed but too weak to hold a bound state, the 
solution evolves from a single unbound pair to a 
very large number of pairs which always have a 
bound state solution. The grand canonical nature 
of the BCS ansatz makes it inherently many-body, 
so a direct connection to the one-pair solution is 
not really possible. 

Five years after the BCS milestone paper, 
RichardsonQ and GaudinQ, succeeded to write 
the exact eigenstate for N pairs interacting via the 
reduced BCS potential, in terms of N parameters 
solution of TV coupled equations. The condensa- 
tion energy obtained from the BCS ansatz, has 
been recovered in the infinite N limit 0, QjJ EH]. 

A decade ago, one of us has developed a new 
framework 12] for many-body effects between com- 
posite bosons. Most of its applications dealt with 
semiconductor excitons. Recently, we have ex- 
tended this framework to Cooper pairs and red- 
erived Richardson-Gaudin equations |l3j. We have 
also succeeded to obtain an analytical solution of 
these equations 14 , 16[ which exactly matches the 
energy obtained through the BCS ansatz. The 
Richardson-Gaudin approach is all the most suit- 
able to investigate the change in condensation en- 
ergy from 1 to N pairs when the pairing potential 
stays constant, because it allows us to handle a 
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fixed number of pairs with Pauli blocking treated 
exactly. 

When one electron pair is added to a system 
already having N pairs, the Pauli exclusion prin- 
ciple shows up in two different ways, (i) Pairing 
(binding) has to use a smaller phase space, so that 
the energy saving per pair due to the attracting 
potential must be smaller in the case of (N + 1) 
pairs than for N pairs since they have less freedom 
to construct the most favorable correlated state. 
This binding decrease is the so-called "moth-eaten 
effect" (ii) In the absence of attraction, the 

additional fermion pair fills the next k-level due to 
Pauli blocking; so its energy also depends on the 
N other pairs. Of course, these two effects must be 
handled self-consistently. However, it is enlighten- 
ing to separate them in order to build some intu- 
ition. The binding energy decrease resulting from 
the "moth-eaten effect" driven by Pauli blocking, 
scales as N/Nq, where N is the number of pairs 
at hands and Nq the maximum number of pairs in 
the potential layer which scales as the sample vol- 
ume. By contrast, the average kinetic energy for 
free pairs scales as e_F, i.e., (iV/L 3 ) 2 / 3 , where L 3 is 
the sample volume in 3D. When N is small, this 
dominates the "moth-eaten effect"; so, the con- 
densation energy per pair, which results from the 
energy difference without and with potential, must 
increases. For larger iV however, the "moth-eaten 
effect" dominates and the condensation energy per 
pair finally decreases. 

This understanding points out an important as- 
pect of the BEC-BCS crossover: It is usually in- 
troduced at the two-body level, a bound-state ap- 
pearing when the attraction passes some thresh- 
old. The many-body solution is then seen as the 
effective potential threshold turning to zero for the 
system to condense at vanishing potential. The 
present work proposes a somewhat different under- 
standing which better bridges 2-body to TV-body 
systems in a BEC-BCS crossover: as the pair num- 
ber increases, a correlated state develops at a lower 
but still finite potential. 

The paper is organized as follow: 

In section [231 we describe the model. We recall 
the single pair case and then turn to a qualitative 
understanding of the many-pair system through 
the condensation energy change when the pair 
number increases. We show that this change is 
quite different in 2D, with a constant density of 
state, and 3D with a density of state which can- 
cels at zero energy. To support this understanding, 
in section [3j we carefully study two pairs through 
the resolution of the corresponding Richardson- 
Gaudin equations with a \fi density of state: we 
show that a binding indeed develops when turning 
from one to two pairs for a potential set exactly 



equal to the threshold value for one pair. We then 
conclude. 



2. Physical understanding 

2.1. The model 

We consider N pairs of fermionic atoms with 
creation operators and ruled by the hamil- 
tonian H = Hq + Vbcs- For same mass atoms, 
the kinetic part Hq reads 



H = X! e k(4 a k 



+ blK) 



(2.1) 



We take as potential the reduced BCS potential of 
standard superconductivity, but without its frozen 
core, namely 
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where = alb_ k while = 1 for < £k < fi 
and zero otherwise; so attraction between zero- 
moment pairs acts from zero to a sharp cutoff 
ft. While this cut-off bears no connection with 
phonon energies, we can still relate it to a physical 
quantity, the scattering length, as shown below. 

2.2. One pair 

The energy E\ of a single pair in this potential 
follows from Cooper's equation 
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(i) In 2D, the density of state is constant. By 
transforming the sum over k into an integral, we 
get for negative E, 
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This function tends to infinity when E — > 0_ and 
to zero as ptt/(—E) when E — > — oo. A bound 
state, solution of Eq. (|2.3I) . thus exists no matter 
how weak v is. It reads E^^ = —2fla/ (1— a) with 
a = er 2 l pv . Note that while p increases linearly 
with sample volume, pv stays constant. 

(ii) In 3D, the density of states can be written 
as p(e) = py/e/Ti where p now is the density of 
state at the potential upper boundary. So 
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tends to p when E — > 0_ and to zero as |pQ/(— E) 
when E — > — oo. A bound state thus exists for v 
larger than a threshold v t h = 1/p. For a poten- 
tial just above threshold, the single pair energy 



tends to zero as e[ 3B) —8(pv — l) 2 Q/n 2 while 
far above threshold e[ 3D ^ rs — |pufi 

(iii) Using this result, we can relate the s-wave 
scattering length a s , commonly used for cold gases, 
to the potential cut-off ft via the density of state 
at this cut-off, p — mL?\/2mVl/2-K 2 . Indeed, for 
fermion pairs interacting via Vbcs, the T- matrix 
for S-wave reduces to 

_ -VUJk . , 

* ~ l-vS(2e k + iO+) [ ' 

wit h, from Eq.(2.5), S" (3D) (2e fe + i0+) ~ p(l + 
vk\J efc/4fi) for 6^ -c ^- The scattering length 
then follows from the scattering amplitude f% = 
—a s /(l + ika s ) which depends on the T-matrix as 
/£ = -mL 3 T°/47r. So, a s ~ mL z v/±n(pv - 1). 
For v slightly above the single pair threshold 1/p, 
we find a s positive with a pair binding energy 
Eb ~ -1/mOj, while below threshold, a s is nega- 
tive and no bound state exists. 

2.3. N pairs 

Richardson 0| and Gaudin Q showed that the 
energy of A fermion pairs interacting via Vbcs 

reads as Em — Ri~\ \-Rn where the iVs follow 

from TV coupled equations 
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(i) 2D systems: Very recently [13, HH, we have 
derived a compact solution of these equations 
when the density of state is constant above a 3D 
frozen core, as for the A Cooper pairs in standard 
BCS superconductivity. Using this result for 2D 
systems which have a constant density of states 
whatever the electron energy is, we get 
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within under-extensive terms in {N/p) n . The 
energy difference without and with potential 
leads to a condensation energy per pair cn = 
[E N {v = 0) - E N ] /A equals 
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the total number of states in the potential layer 
being Nq — J2k Wk = m ^D. This shows 



that 
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decreases linearly with N , due to the 



•'moth-eaten effect" induced by Pauli blocking on 



the number empty states feeling the potential and 
thus available to form the correlated TV-pair state. 
Note that for complete filling, = [ fz^] /p goes 
to zero as 1/p: for N = Nn, the system has lost 
all its freedom to construct a lower energy state. 

(ii) 3D systems: A compact expression of the TV- 
pair energy is not known for a y/e density of states. 
We can yet say that, due to the same decrease of 
available states for pairing, the "moth-eaten ef- 
fect" must bring the condensation energy per pair 
down to zero when N approaches the total number 
of pairs in the potential layer, which in 3D, reads 
TVo = E k ^k = 2pfi/3. 

• For v smaller than Vth{l) = 1/ Pi the condensa- 
tion energy per pair e^ D ' is equal to zero for N = 1 
and also for N = Nq due to the "moth-eaten ef- 



fect". It is however clear that e 
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cannot stay 



equal to zero for all N because when N gets large, 
we can always think of freezing Nq of the TV elec- 
trons as in standard BCS superconductivity. The 
density of states for the other (TVq — TVo) states in 
the potential layer is then finite; So, the (N — No) 
pairs in this layer can always condense, no matter 
how weak v is. For No ^> Nq — No, (a thin shell 
over a large core in momentum space) , the density 
of state in the shell can be taken constant as in 
2D. We can then use Eq. (|2.9|) to estimate the to- 
tal binding energy of these (TV — TVo) pairs which 
only enjoy (TVq — TVo), instead of Nq states, for 
pairing. This gives 
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with a — e 2 ' pv where p is the average density 
of states above the frozen core, p(NoZNq) 1 ^ 3 < 
p < p. By taking for p its lower boundary, we 
find that E {3D) {N ) has a maximum which results 
from a competition when Ao increases, between 
the "moth-eaten effect" and the increase of the 
energy contribution 2<7f2/(l — a) from each avail- 
able states for pairing. This argument shows that 
for v smaller than the threshold for one pair, but 
A large enough, a BCS-like collective state must 
develop with a non-zero condensation energy. As 
a result, the threshold potential vth(N) for the 
appearance of A-pair condensation must decrease 
with A, down to essentially zero when A reaches 
a sizable fraction of the total number of pairs An 
in the potential layer. Consequently, the conden- 
sation energy per pair stays equal to zero 
up to TV = A* reached for v = v th (N*). It 
then increases and ultimately decreases to zero for 
A = TVji (see Fig. [T]), due to the moth-eaten ef- 
fect which always dominates when most of the pair 
states available for pairing are occupied. (Since 
the A dependence of above TV* is not yet 



3 



known, we have preferred not to show it on Fig[TJ) 

• For v above i>t/j(l), a finite condensation en- 
ergy already exists for N = 1. Continuity with 
v < v t h{l) leads us to think that e^ D ^ must hrst 
increase with TV and then decrease due to the same 
moth-eaten effect which dominates condensation 
when N approaches full-filling. To physically un- 
derstand this behavior, we can note that the con- 
densation energy results from a difference between 
the free pair and the correlated pair energy, which 
both increase with N due to Pauli blocking. The 
question then is to understand why, for small N, 
the free pair energy increases faster than the cor- 
related pair energy. When one free pair is added 
at the energy level e, the kinetic energy increases 
by 2/p(e) since the fermion density of state is by 
definition related to the energy increase through 
p(e) — AN/Ae. In 3D, with a y/e density of states, 
the energy change when going from N to TV + 1 
pairs, is thus larger at low energy. By contrast, 
correlated pairs are made with all pair states be- 
tween and SI. When the pair number increases 
from N to N + 1, a very small fraction of each 
of these (0, SI) states is blocked, so that the ki- 
netic energy change is far smaller than when block- 
ing a single low energy state. Consequently, when 
./V is small, the energy difference between adding 
one free pair and one correlated pair must increase 
with N. By contrast, when N gets large, the en- 
ergy cost 1 / p(e) to add one free pair is essentially 
equal to the cost to block most k states making 
the correlated pairs. We are then left with the 
"moth-eaten effect" on the bound state itself and 
the condensation energy decreases. 

We wish to stress that this understanding is fully 
supported by the monotonous decrease found for 
e^ D ^: indeed, the 2D density of state being con- 
stant, the energy cost to block a low energy state 
is the same as the one to block any other states 
making the correlated pairs. The "moth-eaten ef- 
fect" controls the whole behavior and the conden- 
sation energy per pair always decreases when N 
increases. 

We here use a sharp potential cutoff SI for con- 
venience. Real systems do not have such a cutoff; 
so, the fact that the condensation energy vanishes 
when all states below SI are filled, can be ques- 
tioned. A sharp cutoff SI mimics the fall-off of the 
potential at high energy. The potential amplitude 
defines an energy scale, which in turn defines a 
phase space. Beyond this space, the potential is 
too weak to substantially affect the fermion dis- 
tribution which just behaves as free, with a zero 
condensation energy. 

To establish this qualitative understanding on 
stronger grounds, let us concentrate on 2 pairs. 



The change from 1 to 2 pairs will give us the trend 
for the TV-dependence of since, for N = 2, 

Pauli blocking which drives this TV-dependence, is 
already present. 

CM 








1 N* N n N 

Figure 1: Condensation energy per pair as a func- 
tion of pair number: in 2D (dashed line), it always 
decreases due to the "moth-eaten effect" while in 
3D (solid lines for v below or above threshold 
for binding a single pair i>th(l) = this ef- 

fect dominates for large N. As a result, in BEC- 
BCS crossover the condensation energy per pair 
has a maximum which results from the competi- 
tion between moth-eaten effect and kinetic energy 
increase when adding a pair. 



3. Condensation energy for two pairs 

According to Richardson-Gaudin, the exact 
eigen-energy of two Cooper pairs reads as E% — 
Ri + i?2 with (Ri, R2) solution of 

k 

3.1. 2D systems 

We have solved these coupled equations analyti- 
cally for a constant density of states above a frozen 
core [l6| . Using this work for 2D systems with a 
constant density of states p between and SI, we 
get the energy difference between two correlated 
pairs and two single pairs in the large sample limit 
as 

S (2D) _ 2S (2D) = 2 / + 2a_\ + 1 

p\ 1-aJ p 2 

Since 2/p is the exact kinetic energy cost to go 
from 1 to 2 pairs when the density of state is con- 
stant, we find a binding energy decrease when go- 
ing from one to two pairs equal to 2a/p(l — a) in 
agreement with Eq. 12. 81 taken for N — 2, as a 
consequence of the "moth-eaten effect" . 
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3.2. 3D systems 

The situation is more complex in 3D because 
the potential has to be higher than a threshold to 
sustain a bound state for one pair. 

(i) Potential above threshold v > Utfc(l) = l/p. 
A single pair has a bound state with E\ finite neg- 
ative. The energy change between two correlated 
pairs R\ + R2 and two single pairs 2E\ can only 
come from Pauli blocking due to the very peculiar 
form of the reduced BCS potential [3]. We thus 
expect R\ + R 2 ~ 2i?i for L 3 — > 00, i.e., p — > 00. 
This leads us to expand the sum of Eq. (|3.1[) in 
terms of {R x - E x ). Using Eq. (POj) . we get 



= J2(R 1 -E 1 rJ2 



(2e k - E 1 _ 



n+l 



Ri — R2 
(3.3) 

The sum over k, calculated with a py/ e/fl density 
of state, gives 
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where K n (x), defined as 
K n (x) = 



Vydy 



2(y + l)«+i 



-Ex' 
(3.4) 

(3.5) 



goes to a finite value K n when x goes to infinity, 
which is the relevant limit since \E\ \ <C fl. By set- 
ting Ri — Ei{l — ti), Richardson-Gaudin equations 
(OLD then lead to 
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The sum and difference of these two equations give 

00 

= X>x+*2)^« (3-7) 

71=1 

OO 

-A 2 = (ti-t 3 )£(t?-*J)/ir„ (3.8) 

n=l 

where A 2 = (^n)(3j|-) 3 ^ 2 . Two different regimes 
thus appear: 



For A 2 < 1, Eq. 
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Since i 2 + fcj is just 



+ <2) 2 + (*i - t 2 ) 2 } /2, Eq.02J gives (i t 
t 2 ) ~ X 2 K 2 /2K 2 . So the two correlated pair en- 
ergy reads for p large as 
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i?2 thus increases with sample volume as l/p ~ 
1/L 3 . To get the condensation energy change from 



one to two pairs, we must compare this increase 
with the one of two free pairs. Due to Pauli block- 
ing, the second pair momentum must be 2it/L] 
so, the kinetic energy increase for two free pairs 
scales as l/£ 2 , which for large L, is larger than 
the two correlated pair increase. The condensa- 
tion energy per pair, which is the energy difference 
without and with potential, thus increases from 1 
to 2 pairs - in agreement with the continuity argu- 
ment of section 2. 

• The regime A 2 of the order or larger than 1 
corresponds to E\ small, more precisely — E\ < 
/ p 2 / 3 ss 1/mL 2 . The single pair binding en- 
ergy is then smaller than the kinetic energy dif- 
ference between the two lowest free states: in this 
limit, the discrete sum over k cannot be replaced 
by an integral and the above procedure fails. This 
small Ei regime must be handled along the E\ = 
case studied now. 

(ii) Potential at threshold v — l/p. The single 
pair energy E\ then cancels; so we cannot rescale 
the i?i's in terms of E\ as previously done. It 
is reasonable to think that, as for E% small, the 
threshold behavior for E\ — 0, cannot be derived 
by replacing the sum over k by an integral. To 
prove it, let us first subtract the two Richardson- 
Gaudin equations (|3.ip . This gives 
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For R 1 +R 2 = E 2 = 2R and R 1 -R 2 = 2iR' with 
R real but R! a priori complex, we get 
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which also reads 
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This imposes R' 2 real, i.e., (Ri,R 2 ) both real or 
complex conjugate. 

• Let us first show that (Ri,R 2 ) cannot be 
complex conjugate. For that, we use Ea. (|2.3p for 
Ei = 0. Eq.fl3H) then gives 



i?i(i?i — R 2 



Y 2e k (2e k -i? 1 ) 



(3.13) 



If (Ri,R 2 ) were complex conjugate, we could set 
Ri = R 2 = 2Qre ie with r real and < 6 < vr. For 
(Ri,R 2 ) far enough from the real axis to possibly 
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replace the k sum by an integral, we get 
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We expect a 2-pair energy E 2 far smaller than 
the potential threshold, \R\ <C 257, i.e., r <C 
1. The R.H.S. of Eq. (157141 then reduces to 
(m - 2^ cos0/2)(l + 0(y/f)). The real part of 
the above equation then gives — 1 « 4pf2cos 2 0/2 
which is impossible. 

• We now show that (RijR?) cannot be both 
real outside (0,f2). For that, we add the two 
Richardson- Gaudin equations (I3.1[) and use 1/v = 
^ k Wk/2ek at threshold. This gives 



2e k 
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If , R2) were both real outside (0, f2), each term 
of the bracket would be negative; so, the above 
equation is not fulfilled. 

• We are thus left with (R±, R2) both real, with 
one at least in (0, f2) - or possibly (i?i, R2) complex 
conjugate but very close to the real axis so that 
the k sum cannot be replaced by an integral. To 
estimate the two pairs ground state, we can reduce 
this discrete k sum to its first few terms, namely 
k = 0, i.e., 6k = and (k^k^kj) = ±2ir/L, 
i.e., £k = (±27t/L) 2 /2to = ej,. The Richardson- 
Gaudin equations with these seven terms in the 
sum reads at threshold 
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The proper rescaling then is Ri — 2e£it;. So the 
two Richardson-Gaudin equations read 
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where peL scales as L. 

- We first look for (1*1,112) both real. For L 
large, we get —1/ui ~ 2ptL ~ 6/(1 — u 2 ). This 
gives i?i = —1/p and R 2 = 2el — 6/p; so, the 
two-pair energy reduces to R\ + R 2 ss 2e^ — 7/p . 

- We then look for (ui,u 2 ) complex conjugate, 
i.e., ui = u + iu' = u\. It is possible to show 



that Eas. (|3.17p then have no solution for u < u' 
or u ~ v! while a solution exists for v! > u. It 
reads u w 1 — 5/2peL with u' cc l/pej,. This leads 
to a 2-pair energy i?i + i? 2 ~ 4e^ — 10/p which is 
above the energy found for (Ri,R 2 ) real. 



Since the kinetic energy of two free pairs is 2cl , 
the condensation energy for two pairs at the po- 
tential threshold for one pair , v — 1/p, reduces to 
2e£ — (2e£ — 7/p) = +7 / p. Since this condensation 
energy is zero for one pair, it indeed increases from 
N = 1 to ./V = 2 in agreement with the continuity 
argument of section 2. As a result, the potential 
^tft(l) = 1/Pi too weak to sustain a bound state 
for a single pair, can bind two pairs: by continuity, 
the potential threshold for pairing, v t h(N), must 
decrease when the pair number increases. 

4. Conclusion 

We show that a system of fermion pairs with 
attraction too weak to sustain a bound state, 
can nevertheless form a many-body correlated 
state when the pair number increases. Us- 
ing Richardson-Gaudin equations, we demonstrate 
that the required potential strength for conden- 
sation decreases as the pair number increases. 
For strong attraction, two fermions form a bound 
state; the many-body state can then be viewed 
as a congregation of single pairs with Pauli block- 
ing decreasing their binding energy ("moth-eaten 
effect"). However, when attraction gets weaker, 
more than one pair is necessary to form a con- 
densed state as explicitly shown by going from 
one pair to two pairs at the potential threshold for 
binding a single pair. For a vanishing interaction, 
a large number of pairs are required to obtain BCS 
condensation. A cross-over then appears when the 
pair number increases at constant but finite poten- 
tial. 

This work also reveals, through difference be- 
tween 2D and 3D behaviors, the subtle interplay 
which exists between Pauli blocking at the free 
fermion level which increases the kinetic energy 
when the pair number increases, and Pauli block- 
ing at the paired level through the "moth-eaten 
effect" resulting from the decrease of the number 
of empty states available for pairing. 
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